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Abstract 

We consider classical particles coupled to the quantized electromagnetic field in the background of 
. a spatially flat Robertson- Walker universe. We find that these particles typically undergo Brownian 

00 ! 

, motion and acquire a non-zero mean squared velocity which depends upon the scale factor of the 



universe. This Brownian motion can be interpreted as due to non-cancellation of anti-correlated 
vacuum fluctuations in the time dependent background space-time. We consider several types of 
coupling to the electromagnetic fleld, including particles with net electric charge, a magnetic dipole 
moment, and electric polarizability. We also investigate several different model scale factors. 
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I. INTRODUCTION 



Brownian motion of a particle in a thermal bath is a well-known phenomenon. (See, for 
example, Ref. [l|.) In this case, the particle's mean squared velocity grows linearly in time 
until dissipation effects become important, after which it approaches a non-zero equilibrium 
value. The linear growth phase is characteristic of any random walk process, in which 
each fluctuation is independent of previous fluctuations. Quantum fluctuations are quite 
different from thermal ones in that the former are strongly correlated or anti-correlated. 
This does not, however, prevent quantum Brownian motion, which will be the topic of this 
paper. The existence of Brownian motion in the Minkowski vacuum state is controversial. 
Although conventional quantum electrodynamics suggests that the only effect will be an 
unobservable mass renormalization, Gour and Sriramkumar have argued that there could 
be an observable effect on charged particles coupled to the fluctuating electromagnetic field. 
Brownian motion in the presence of boundaries is less controversial, and has been studied 
by several authors [1, II, [5j_l, 0, Barton [3.] was the first to examine fluctuations of the 
Casimir force. Wu et al ^calculated the IBrownian motion of an atom near a perfectly 
reflecting plate due to fluctuations in the retarded van der Waals force. The mean force here 
is the Casimir- Polder force [9]. The analogous Brownian motion of a charged particle near 
a reflecting plate was treated by Yu and Ford [^. In all of these cases, the mean squared 
velocity of the particle approaches a constant even in the absence of dissipation. This is 
required by energy conservation, as there is no energy source in these static configurations. 
The fact that the late time mean squared velocity is non-zero can be attributed to the effects 
of switching when the interaction is turned on. Switching effects were recently discussed by 
Seriu and Wu [8j. 

The mechanism which enforces the lack of growth of the mean squared velocity can be 
understood as anti-correlated fluctuations. A charged or polarizable particle in a Casimir 
vacuum can acquire an energy E from a fluctuation. However, that energy is typically 
surrendered on a time scale of order h/E to an anti-correlated fluctuation. The correlation 
functions of the quantized electromagnetic field automatically enforces the required anti- 
correlations [10]. The quantum fluctuations of the stress tensor in flat spacetime also exhibits 
subtle correlations and anti-correlations, as is discussed in Ref. jTl| |. 

The Brownian motion of test particles is an operational means to describe a fluctuating 
quantum field. This approach can be used to treat the quantum fluctuations of the gravita- 



tional fieldjwhich has been a topic of much interest in recent years [12, ll3|, [IJ, Il5|, [Ig, ll7|, 
[i8|,[Ti,[2a,|2l|,[2l,[2l,[2l,[2|. 

In the present paper, we will investigate the Brownian motion of various types of particles 
coupled to the quantized electromagnetic field in the background of a Robertson- Walker 
spacetime. Here the time-dependent background geometry can act as an energy source, so 
the particles can acquire a net kinetic energy. 

In Sect. [Til we develop the basic Langevin equation formalism for calculating the mean 
squared velocity of classical particles coupled to a fluctuating force in a spatially flat 
Robert son- Walker background. The for mal ism is applied to a several specific choices for 
the scale factor of the universe in Sect. IIIII Our results are summarized and discussed in 
Sect. [IE 

Unless otherwise noted, we work in Lorentz-Heaviside units with h = c= 1. 



II. BASIC FORMALISM 

The equation of motion of a classical point particle moving in a curved spacetime with a 
four-force is 



r^™-^. (1) 
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where m'^ is the 4-velocity of the particle, m is its mass and r the proper time. The operator 
D/dr on the right-hand side of Eq. ([1]) is the covariant derivative given by 

Here we take the space-time geometry to be that of a spatially flat Robertson- Walker 
universe with metric 

ds^ = -dt^ + a^{t){dx^ + dy'^ + dz^) , (3) 

where a{t) is the scale factor. We will restrict our attention to the case where the particles 
are moving slowly with respect to these coordinates, in which case the particle's proper time 
becomes the coordinate time t. Because of spatial isotropy, we can consider a particular 
direction, the x-direction, and write the equation of motion as 

^ + 2V = lr. (4) 

dt a m 

Here we have used V^^ = T^^ = d/a, where d = da/dt. We will take the four-force to be of 
the form 

r = /'^ + , (5) 

where f^^^ is a non-fiuctuating external force, and f'^ is a fluctuating force produced by the 
electromagnetic vacuum fluctuations whose mean value vanishes: 

{f'n = . (6) 

First, let us consider the case of free particles, which corresponds to the case when 
fext = 0. Thus Eq. (jlj) can be written as, 

1 d f'^ 
dt ^ m 

which after integration reduces to 



'itf)u-{tf) - a\to)u-{to) = - / dta'{t)f'-{t). (8) 



m 



to 



Assuming that these particles are initially at rest {u^{to) = 0), we find that the velocity- 
velocity correlation function is given by 

{u%tf,n)u^itf,r2)) = [ dt,dt2a\t,)a\t^) {f''{t,,r{)f'%t2,r^)) . (9) 

TDj d J J 

Another case of interest is when there is an external force which cancels the effect of the 
cosmological expansion: 

/:.. = 2-«^ (10) 

a 

This is the case for any particles in bound systems such as galaxies or molecules. Such 
particles do not participate in the cosmological expansion and in this case two such particles 
do not move apart on the average. We will refer to these as bound particles. In this case. 
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and the velocity correlation functions for particles which start at rest at t = t/ is 

{u%tf, n)u%tf, r,)) = —! dhdh (/'^(ti, n)f'%h, r^)) . (12) 

Tfl J 

Note that the above expression is a coordinate velocity correlation function. In a 
Robert son- Walker space-time, proper distance between particles, If, is related to the co- 
ordinate separation r at t = t/ by // = a/r. Thus the proper velocity correlation function is 
given by 

{v-{tf,n)v'\tj,r,)) = a}{u%tj,n)u-{tf,r,)) . (13) 



A. Charged Particles 

In this section, we will consider electrically charged particles with charge q coupled to a 
fluctuating electromagnetic field. In this case, the four-force is 

fx ^ Q_pxt^^^ _Q_pxt _ (14) 
m m 

For the case of free particles, Eq. yields 

{u%tf,n)u%tj,r2)) = ^Jdt,J dh a\h)a\t2) {{F^'ut), {F^-'ut)^)Rw, (15) 

where the sub-indexes 1 and 2 refer to the coordinates (ti,ri) and {t2,r2), respectively, and 
the subscript RW denotes a vacuum correlation function in Robertson- Walker spacetime. 

This correlation function is obtained from the corresponding correlation function in fiat 
spacetime by a conformal transformation. First write the Robertson- Walker metric in its 
conformal form 

ds^ = a^{-drf + dx^ + dy'^ + dz^), (16) 
with dt = adrj. The field strength tensor in these coordinates is given by 

(F^Rw = (F^M , (17) 

where the subscript M refers to the Minkowski space fi eld strength. This may be seen, 
for example, from the fact that the Lagrangian density, yJ—gF^^'^F^j^ is invariant under the 
conformal transformation. From this and Eq. f[T^ . we find 

(u^(t;,n)n^(t/,r2)) = -^Jdmj dm{F^%m,r{)F^^{m,r2))M. (18) 

Here the appropriate component of the Minkowski space correlation function is given by 
Eq. (IA5I) . The key feature of this result is that the scale factor does not appear inside the 
integrand. Thus the cosmological expansion has no effect on the Brownian motion, and 

hence we do not find an interesting result in this case. 

The case of bound charged particles is different. In this case, from Eqs. (fT^ and (lASp . 
we find 

{u%r,,n)u%7^,r2)) = ^Jdr],J dr,2{{F'^ Ut) ^ {F''ut)^)nw (19) 

Here r is the spatial separation of the particles and r'^ = — 2Ax'^. Now there are factors 
of in the integrand, which will lead to non-trivial effects in an expanding universe. 
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B. Magnetic Dipoles 

In this section we will consider particles with a magnetic dipole moment. In flat spacetime, 
such particles experience a force when there is a non-zero magnetic field gradient: 

/ = -Vm (20) 

where u = —fl.B, is the magnetic potential energy, fi is the magnetic moment and B is the 
magnetic field. Writing this force in covariant form, we have for the x-component 

r = fid' {F^'x^x'^) (21) 

where, = (0, 0, 0, a) and x'^ = (0, 0, a, 0), and such that x^x^ = x'^x'^ = 1. 
For free particles, the velocity- velocity correlation function, Eq. (Q, becomes 

(M^(t;,ri)n^(t/,r2)) = ^i^^^ J dt, j dtA.d,, {{F^'^z^'y), {F^'x^x'^) ^)nw. (22) 



Again, we may use Eq. f[T7j) to write the above expression in terms of 
dx^dx2{F^'^ {rji^rij F^y{r]2,r2))M, which may be evaluated to write 

(M^(r/,ri)M^(77,r2)) = - ^2q4(^^) / / dV2a-\rii)a'\m) x (23) 



^ 6(r/2 - Vi 



\2 



In Eq. ( l23l) . we have used the coincident limit r — > only in the numerator, in order to 
simplify our expressions. This will not alter our final results, because we will take this limit 
after the integrations. 

For bound magnetic dipoles, we may start with Eq. f[T^ and follow the same procedure 
to find 

(M^(77,ri)M^(77,r2)) = [ drii [ dri2a~^{rii)a-^{r]2) X (24) 



2 6(?72 - rji) 



2 



+ 



where again the coincidence limit in the spatial coordinate was taken in the numerator of 
the integrand. 



C. Polarizable Particle 



We will consider in this section a polarizable particle, described as a point particle with 
a static polarizability a. In an inhomogeneous electric field, such a particle experiences a 
force 

fix) = |V(E^) , (25) 

which in covariant notation becomes, 

r = ^d"" {F'^UtXxY , (26) 
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where the low velocity limit was taken with u^, = —5°, = (0,a,0,0). For a free particle, 
Eqs. © and lead to 

(w^(t;,r0n^(t/,r2)) = d^, j dr^, x (27) 

where we used the fact that = —aS"^ and = (0, a, 0, 0), in conformal coordinates. 

Here we use the Wick theorem to calculate the two-point function: {E'^{ri)E'^{r2)) m = 
{E,{n)E\r2)Ej{r^)E^{r2))M, finding 

{E,{r{)E\r2)E,{r{)E^{r2))M = 2 [{Ei{r{)E,{r2)) M{E\n)E^ {r^)) m] , (28) 
or, from the procedure outlined in Appendix A, 

{EM)E\r,)E,in)E^ir,)),, = { ,2[!|'^^ I '^^^ ^ ^^p • (29) 

To simplify our expression we will consider the coincident limit in the spatial coordinate 
r = r2 — ri = 0, only in the numerator of all factors, as in the previous cases. Then the 
velocity-velocity correlation function is. 



(u^(77,ri)M^(r/,r2)) = J dr/i J d-q^a ^{j]i)a ^(r^s) x (30) 

32 , 136(r72 - r7i)2 144(r/2 - r/i)^ 



+ 



III. SPECIFIC UNIVERSE MODELS 

In this section we will apply the basic formulas obtained in Sect. |TT] to investigate the 
influence of different scale factors on the Brownian motion of particles induced by quantum 
vacuum fluctuations of the electromagnetic field. 



A. Asymptotically Static Bouncing Universe 

The study of bouncing universe was considered by some authors in the past [2^ . Here 
we will study a special case, which is asymptotically static in the past and future. We will 
take the scale factor to have the form 

= ^2 + C'4" ^^^^ 

where G and ?7o are constants, and n is a positive integer. Note that when r] —* ±oo, the 
universe is asymptotically fiat and a goes to unity. It will be convenient to consider different 
choices of n for different types of particles in order to simplify the corresponding integrals. 
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1. Bound Charged Particles 

In this case, we set n = 1, so that Eq. f|T9l) becomes 



2 r+oo f+oo / 2 , /y2 2\2 

{/\u{^, r,)Au{v, T,)) = ^ dr^J d^, (32) 

J~oo J -oo \ Vl Vo / 



where we used Eq. ( IT9|) . This integr al is evaluated in Appendix B, with the result being 
a rather complicated expression, Eq. flB4p . This result simplifies considerably in the limit 
that ri — > r2 and G ^ 1 to 

(A..=> = (33) 



Note that because aj = 1, this expression also gives the mean squared proper velocity (At>^). 

We can gain some insight into this result by writing it in terms of a characteristic measure 
of the maximum curvature. Consider the scalar curvature R and evaluate it using the scale 
factor Eq. (l3Ti) with n = 1. We find the following result, 

i?o = ^ (1 - G^) , (34) 
'/o 

where Rq is the Ricci scalar when rj = 0. If we consider the limit G ^ 1, i?o is negative and 

Rl « (35) 

'/o 

Then (Af ^) in terms of Rq is 

(At;2) = Z^!i^. (36) 
\ ' 6144m2 ^ ^ 

We can also write Eq. in terms of the redshift, defined by: a"^ = (1 + z)^, with 
= a(0)/a(oo) = (l/G)^, where is the minimum scale factor. Considering 2; ^ 1, we 

get 

(At;2) = Ifl^z^. (37) 
\ / 6144m2 ^ ^ 

The mean squared velocity is proportional both to the squared redshift and to the maximum 
curvature. This can be associated with an effective temperature using the non-relativistic 
equation: ksTeff = m{Av'^). Where Tg// is the effective temperature and ks is the Boltz- 
mann constant in Lorentz-Heaviside units h = c = 1. Then, T^jf is, 

T.n^'-^(¥)'^^ (38) 



where /c = 1/ a/|-Ro| is the lenght curvature and Ac = l/m is the particle's Compton wave- 
length. 
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2. Free Magnetic Dipoles 

Again take the scale factor to be Eq. fl3Tl) with n = 1. Then Eq. fl23|l for the mean 
squared velocity becomes, 



-2 Q{m-r,if 



Using the same procedure as before, we find the following result in the coincident limit when 
G'> 1 



,2 

In terms of the redshift (2; S> 1), we have 



(A.^) = l^M. (41) 



In contrast to the result for bound charges, the effect decreases with z. For the case of 
electrons, it is convenient to write the magnetic moment as 

/i ~ ^ . (42) 
The effective temperature in terms of curvature length and Compton wavelength Ac is. 



3. Bound Magnetic Dipoles 

Here we choose the scale factor to be of the form of Eq. fl3Tl) with n = 3. In this case, 
the mean squared velocity from Eq. (|2^ is. 



^2 r-oo roo / 

{Au{ri,r])Au{r2,r])) = — / dr]i dr]2 i 

m j_oo -'-00 V 



6(?72 - ?7i)^ 



(44) 



We may evaluate this integral using the same technique as before, with the result in the 
coincident limit 

{Au^) = ^V^^ r^. 45 

The physical velocity when G ^ H is, 

) = -^i^^^l^ol ^0, (46) 
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where now the scalar curvature at = is 

Ro = --^G'l^. (47) 
In terms of the redshift, given by 1 + z ^ 2 = a(oo)/a(0) = G^/^, we can write 

which shows that the effect of quantum fluctuations grows with z. Associating an effective 
temperature we have, 

10- V /A 



Here the temperature grows with z as in the bounded electric particle due the extra force 
that acts on the magnetic dipole. Indeed, we see that the effect here is smaller than that 
one indicated by Eq. ( l38l) . 



4- Free Polarizable Particle 

Again we take the scale factor to be of the form of Eq. (|3T1) with n = 3. Equation (|3T|) 
for the mean squared velocity can be written as 

(AM(ri, r/) AM(r2, r/)) = ^^,^4^, x (50) 

J-oo V vl + vi J \ vi + vi J 

32 136(r/2 - Vif 144(7^2 - Vi] 

r / \ o . m c 



Following the procedure previously used, we find, in the coincident limit, 

2^ 3X 10-2^2 - 1)' 

This physical (comoving) velocity can be expressed in terms of Rq given in Eq. fHTj) as 

At; = 2 2 — ^. (52) 
We can also write Eq. (1521) in terms of the redshift as 

The mean squared velocity decreases with the redshift in contrast with the bounded 
particle cases investigated in the previous section. This is due to the fact that the atoms 
are free of external forces. This effect can be associated with an effective temperature using 
the non-relativistic equation: fc^Te// = m{Av'^). Thus, we obtain 

10-V fxX 2 

This result shows that the temperature decreases with the redshift because the particles are 
free of external forces. 
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B. Asymptotically Bounded Expansion 

A universe with asymptotically bounded expansion was studied in Ref. ^3] 5 where the 
production of massive particles were considered. Here we will investigate the Brownian 
motion effects in scale factors of the form 

a" = Oq + a" tanh(?7/?7o), (55) 

where n is a positive integer, ao and ai are dimensionless constants and r^o is a constant with 
dimension of time. We note that when t] — >■ ±00 ^ — > Oq ± a^. Then, this universe is 
asymptotically flat in past and future, but it is not symmetric and exhibits only expansion. 



1. Bound Charged Particles 

Here we take the scale factor to be given by = 2 in Eq. (l55l) . The mean squared 
coordinate velocity is then given by 



^2 r+co 

{Au{ri,ri)Au{ri,r2)) = ^ dr]2 



drji 



X 



Oq + af tanh 



al + af tanh (^^ ^ 



(56) 



This integral is calculated in Appendix B, with the result given by Eq. ( IBSp . 

In this model, the physical velocity is related to the coordinate velocity by (Af ^) = 
aj{Au'^), where aj = Og + . If we will make a Taylor expansion in r up to the zeroth order 
term, we find 



-4g2(ag + aj) sinh" 



9-^ + 3C(3)) , 



(57) 



where a = ao/ai. Here the expression in parenthesis is a negative constant and ({x) is the 
Riemann zeta function. We can write Eq. (157|) in terms of the scalar curvature at = 0, 
given by 



Rn 



(58) 



It is also interesting write the mean squared velocity in terms of the redshift defined here as 
1 + z ^ z = a{oo)/a{—oo). Then, (Af^) in terms of Rq and the redshift is given by 



{Av') ^ 10''^Roz\ 



(59) 



when z ^ 1. The effective temperature is now 



eff 



10 



-2 r 



z\ 



(60) 
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2. Free Magnetic Dipole 

Now take the scale factor to be Eq. fl55l) with n = 1. Then the mean squared velocity is 
given by 



{Au{ri,r])Au{r2,r])) 



X 



1 



X 



+00 r+oo 

oo J^oo ^^"^ V«o + fli tanh(r7i/r7o) J \ao + ai t&nh{r]2 / r]o] 
-2 6(?72-^?i)^ 



Following the same procedure as before we find that in the coincidence limit (r 



{An') 



sinh 



2 Va - 1 



[C(5) - C(6)] 



As before, a = ao/ai. The physical velocity is 



where now 



Rn 



6al 



The temperature in terms of Ac and Ic is, 



Teff ^ 10 



^3 r 



(61) 



0), 



(62) 



(63) 
(64) 

(65) 



3. Bound Magnetic Dipole 

In this case, let n = 3 in Eq. fl55p . The mean squared velocity becomes 



{Au{ri,ri)Au{r2,T])) 



X 

1 



X 



-foo p+oo 

oo '^^ J-oo V^o + aitanh(?7i/?7o)y V^o + cti tanh(?72/?7o, 
-2 6{V2-Viy 



[-(^?2 - ^71)^ + T^]^ [-{r]2 - r]i)'^ + r"^ 
In the coincidence limit (r ^ 0), 



2, 3 X 10-V^o^^ 



[C(5)-C(6)], 



(66) 



(67) 



where 
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The effective temperature is given by 



4-- Free Polarizahle Particle 

The scale factor is again given by Eq. fl55l) with n = 3 and the mean squared coordinate 
velocity is given by 



(Au(r7,ri)Au(r7,r2)) 
1 

a3 + a?tanh(2i 



a 



dr]2 



drji X 




32 



a3 + a?tanh(^2i^ 
136(r/2 - r/i)2 



144(772 - //I; 



I [-{V2 - ViV + r^V [-iV2 - ViY + r^f ["(^2 - Vif + r^V 
Following the method used previously, we find the result 



(70) 



{An') 



lOa^ 



sinh 



In 



a 



a-' 



X [C(9) - C(10)] 



(71) 



in the limit r — >^ 0. Here C,{x) is again the zeta function, and the expression in brackets is 
positive. The physical speed satisfies the relation 



lOa^ 



sinh 



2 



(72) 



We can write Eq. (172|) in terms of the scalar curvature Rq = 2a1/{r]Qal), and of the redshift 

as 



{Av^) ~ 10 



-11 d4^4 



(73) 



when z ^ 1. The corresponding effective temperature is 



eff 



10 



-11 



a 



A. 



(74) 



Note that for this class of scale factors, we find the unexpected result that (Af ^) oc for 
all fo ur t ypes of particles being considered. These results will be discussed in more detail in 

Sect, nvi 
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C. Another Bouncing Universe 

Here we will consider universes with scale factors of the form 

a- = H\rf + r,l), (75) 

where n is an integer, if is a constant with dimension of inverse of time, and rjQ is also a 
constant but with dimension of time. Although these models are asymptotically fiat in the 
past and in the future, the scale factor does not approach a constant, in contrast to the 
models in Sect. IIII A[ 



1. Bound Charged Particles 

Here we take n = 2 in Eq. (175]) . In this case, the mean squared velocity expression is, 

(An(,, rO A.(,, r,)) = - dr,, dr,, ^^^^^^^^ J (76) 



1 

H'^r]2 + H'^i]q J 7r2[— (r/2 — r7i)2 + r^]^ 

This integral is evaluated in the Appendix, resulting in Eq. (1B13I) . If we take the r limit 
of this expression and use the fact that aX rj = rjf ^ rjQ, the scale factor is a{rjf) ^ Hrjf, we 
find 



.2, _ . 



The Ricci scalar curvature when rj = is 



6 



Ro = ■ (78) 



In terms of this curvature and the redshift, 1 + z ^ z = a//a(0) = rif/rjo, when rjf ^ rjo the 
mean squared velocity turns into 



Teff-^^il^] z\ (80) 



The effective temperature is. 



Comparing Eq. (IHOl) with Eq. (!38l) we see that they are the same except for a numerical 
factor. 
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2. Free Magnetic Particle 

Consider now the scale factor given by setting n = 1 in Eq. fl75l) . We find that the mean 
squared velocity is, 

(Au(ri, r/)AM(r2, t])) = x (81) 



r"00 

drix I dri2 



oo 



-2 Q{r]2-ilif 



[-{ri2- riiY + r'^Y [-(^2 - ^71)^ + r^] 
In the coincidence limit 

and the physical velocity is, 

Af = z , (83) 

with Rq = — 12/(/7^?7o). This is essentially the same result as in Eq. (|4T1) . 



5. Bound Magnetic Dipoles 

Consider now the scale factor obtained by setting n = 3 in Eq. (!75l) . In this case, the 
mean squared velocity is 

{Au{ri,r])Au{r2,r])) = — X (84) 



00 poo 



00 J —00 



2 6(772 -r?i)^ 



-("'72 - ^7i)^ + ^^]^ [-('72 - ''?i)^ + 
In the coincidence limit 



, ^ o, 6 X 10 "^u? , ^. 



with the physical velocity now given by 



where here 



/ A 2\ 3 X 10 V^|i?oP 2 /Q«^ 



and 2; = a(r7j)/a(0) = {f] j / rjo)'^^^ ■ This is of the same form as in Eq. 
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4- Free Polarizable Particles 

Now let us take the same scale factor as in the previous subsection, Eq. fl75l) with n = 3. 
In this case, the velocity-velocity correlation function is 

+ 00 /•+00 -1 -1 



oo 



32 ^ 136(7^2 -r7i)=^ ^ 144(r/2-77i: 



Doing the two integrations and taking the limit r — * 0, the mean squared coordinate velocity 
is 

, ^ 2^ 3 X 10-2^2 

and the proper velocity is nearly the same as in Eq.f l52l) . 

3 X 10-2a2 

where now 



(Av^) = \ ^ \Ro\'z-\ (90) 



D. Oscillatory Expansion 

In this section we add an oscillatory term in the bouncing electric particle case studied in 
the last section. This is an oscillatory expansion universe which was treated in some works 
as, for example in [28|. Here we will see the effects of the amplitude of the oscillation in the 
mean squared velocity for the charged particle case. The scale factor is taken to be 



a 



2 



4iv'' + vl) + alcos{ujri), (92) 



where oq and the frequency u are constants with dimension of the inverse of time, the 
amplitude ai is a dimensionless constant, and rjo is a constant with dime nsion of time. Thus 
the oscillatory term is a small perturbation of the case studied in Sect. IIII C II The mean 
squared velocity is 



q2 r+oo p+co 

{Au{ri,ri)Au{ri,r2)) = — dr]2 j drj 



-{V2 - Vif - r 



2 



X 



1 \ / 1 



(93) 



alr]f + alr]Q + aj cos{ujr]i) J V'^o'?! + (^oVo + '^i cos{ur]2 
In Appendix B this integral is evaluated, and can be shown to lead to the result 

{Av{ri)Av{r2)) = {Av )o H {3uriosmh{ujr]o) + cosh(u;r7o)) , (94) 
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where // = ajr, and {Av^)o is the oi = result found in Eg. (1791) . Notice that the velocity- 
velocity correlation function in this model will diverge in the coincidence limit. So, we 
can not obtain the mean squared velocity in this limit. This reflects a breakdown of our 
model in which the particles are treated as classical point objects. Our model requires that 
// ^ Ac, where Ac is the electron Compton wavelength. In any case, our perturbative result 
requires that the second term in Eq. flMl) be small compared to (Af^)o. Nonetheless, we 
can conclude that the oscillations tend to increase the mean squared velocity in a way that 
grow exponentially with u in the limit that utjq ^ 1. 



E. de Sitter Space 

In this section, we will investigate the effects of the vacuum fluctuations in de Sitter 
space-time. It is well known that de Sitter space can be considered as some special stage of 
the universe history, which is known as the inflationary phase of the universe. It was Guth 
(29l | who first noticed that using some exponential expansion of the universe it would be 
possible to solve three of the standard universe's modelproblems: 1) the flatness problem, 2) 
the horizon problem and 3) the primordial monopoleproblem. This scenario was extensively 
developed since the Guth's original work (see some good reviews about inflation in jSOj], [3ll |) 
and n owadays it seems to be in good agreement with the observations [32|, [33l |. 

We use a scale factor in the form: 

^ (95) 



where — cx) < ?7 < 0. We restrict our attention to the range rji < rj < rjf, where |?7/| <^ |?7i|. 



1. Bound Charged Particle 

Here the velocity-velocity correlation function is given by, 

{Au{r],ri)Au{ri,r2)) = ^ dr]i dr]2{HT]if{Hr]2 

_7r2[-(?^2 - ^lY + r'^Y _ 



(96) 



Notice that the upper value in the integral range is put to be zero to simplify our calculations, 
but in fact it is very small but not null. Using Maple, we find the following result. 



{Au{r,)Au{r2)) = ^ [H ' 



If \rii\ ^ r the coordinate velocity is: 

{Au{ri)Au{r2)) = 
and consequently the physical velocity is, 

{Av{ri)Av{r2)) 



9 9 * 



iv?l. 



2/2 • 



(97) 



(98) 



(99) 
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Our answer is positive and constant and as in the oscillatory case it depends on the initial 
proper particle separation, /j. If k ^ Ac, then Eq. (p9|) reduces to. 



{Av{ri)Av{r2)) - j - 10"^. (100) 



2. Free Magnetic Dipoles 

Using fl95l) and fl23l) we obtain the following mean squared velocity, 

(AM(ri,r/)AM(r2,r/)) = ^f^^ / rfr/i / dr]2 iHr]2) (Hiji) x (101) 

2 



4 2 



-2 6(?72 - r^i) 
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which, after integrations results in. 



(An(ri, r^) An(r2, r^)) = , (102) 

From the equation above it is not possible to get the coincidence limit because the r^- 
divergence. But we know that If = ajr, then, 

{Avin, r])Avir2, v)) = (103) 

Considering q/m, m 1/Xc and If Xc, we have 

{Av{ri, 7])Av{r2, v)) ^ ^ (104) 

The effect in de Sitter universe does not depend of the time. We could also make a estimate 
of the constant H , which has an inverse of length dimension . If ~ 1/Z, we have, 

2 j)^2 

{Av{ri,7])Av{r2, v)) - (105) 



Assuming that / ~ Ac, we get 



g2 

(At;(ri,r7)At;(r2,?7)) ~ — , (106) 



which is basically the same result we found for the bound charge case. 



3. Bound Magnetic Dipoles 

Using the scale factor given by Eq. (1951) and Eq. (1241) . we find that the mean squared 
velocity is given by 

,2 rO /"O 



{Au{ri,r])Au{r2,v)) = ! di^i ! dr]2 {H7]2f {Hr]if x (107) 

2 , Q{V2-Vi? 
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In the hmit which r is very small, we have 



(AM(ri,r/)AM(r2,r7)) = "4^. (108) 



Thus, the physical velocity-velocity correlation function is given by 



{Av{n, r])Avir2, v)) = (109) 



which is negative. Making the same estimates as in the free magnetic particle case for H, 
we obtain, 

{Avin,r])Av{r2, v)) = (HO) 

Negative mean squared velocities have been found by previous authors P, H, [1] , and can 
be interpreted as a reduction in quantum uncertainty. It is well known that a quantum 
massive particle is described by a wave packet which must have a position and momentum 
uncertainty given by the Heisenberg uncertainty principle (Ap^Ax) > 1. If the uncertainty 
in position is such that Ax < /, where / is the average separation between two dipoles, we 
have, {Avcc)q > l/{lm), or 

(AO,>f, (111) 
which is larger than the magnitude of the right-hand side of Eq. flllOp . 



F. Radiation Dominated Universe 



The radiation dominated era of our universe is usually defined as the early period when 
radiation and relativistic particles were usually more important than ordinary matter. Here 
we will investigate the effects of this important universe stage in the average squared velocity 
and evaluate the effective particle's temperature in the beginning of this era in the bounded 
electric particle case. Consider the scale factor. 



Then the mean squared velocity is given by 



{Au{r],ri)Au{r],r2)) 



no 



-iV2-Vi) 



^12 



Taking into account the condition r <^ r/o; this integral results in 



{An') 



where Ci ^ 6 x 10 ^ . The physical velocity can be written as 



(112) 



(113) 



(114) 



;ii5) 
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The redshift factor here is 1 + z ^ z = rjf/rjo. The scalar curvature vanishes for this metric, 
but a reasonable measure of the characteristic curvature is a typical component of the Ricci 
tensor in an orthonormal frame, which gives 

^0^77^. (116) 



Thus, we can write 



{Av') = ^Roz\ (117) 
which is essentially the same as the result found in previous cases in Eqs. f l371) and (I79|) . 



G. Matter Dominated Universe 



Now consider bound charged particles and a scale factor of the form, 

a' = Hy. 
The mean squared velocity in this case is 



{Au{r],ri)Au{r],r2)) 



2 /-oo 



dr]i / dr]2 



1 



1 



-{r]2 - Vif - r'^ 

Performing the integrals in the limit r ^ r^o, we obtain the following result 



2, 10-5g2 



The physical velocity is given by 



/A 2\ C2g' ^/ C2q' „ 2 



;ii8) 



(119) 



(120) 



;i2i) 



where C2 — 10 ^, z = {rif/rji)'^, and Rq = 1/(11*1]^). Again this is of the same form as 
Eq. dSZD. 



IV. SUMMARY AND DISCUSSION 



We have investigated the Brownian motion of particles coupled to the electromagnetic 
vacuum fluctuation in Robert son- Walker universes. We considered several types of particles, 
including ones with electric charge, a magnetic dipole moment, and electric polarizability. 
We also allowed both the possibility that the particles are free, moving apart on the average 
as the universe expands, or bound by a force which cancels the effect of the expansion. 
Our results for the mean squared velocity induced by quantum fluctuations can typically 
be written in terms of a characteristic measure of the space-time curvature, Rq, and a 
redshift factor z. Our treatment assumes semiclassical point particles which always move 
non-relativistically in the comoving frame. Thus our results are restricted to cases where 
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(Av^) -C 1. We should also note that we are working in a regime where quantum particle 
creation by the gravitational field [35j] is small, so we should also requir e that Ro /m? -C 1. 

In many cases, such as the classes of scale factors studied in Sects. IIII Al and IIII Ct the 
effect for bound particles tends to grow with increasing redshift factor, where as that for free 
particles goes the other way. This might be due to the fact that the bound particles are not 
on the average moving on geodesies, and be subject to an acceleration radiation effect of the 
type first studied by Unruh |3J]. However, in other case , such as those treated in Sect. IIIIBl 
the effect scales differently with redshift. These results need to be better understood. 

In most of the cases studied, the m ean squared velocity is finite and positive. However, 
in a few cases, such as in Sects. IIIIDI and IIII E 31 , we found a velocity- velocity correlation 
function which is singular at spatially coincident points and can be negative. Both of these 
phenomena signal a breakdown of our approximation of point classical particles coupled to 
a quantized field. 

However, our view is that these results still have physical content if properly interpreted. 
The spatial separation of particles should always be large compared to the Compton wave- 
length, and the separation should be sufficiently large to insure that (AfiAf2) <^ 1. With 
these restrictions, one can still conclude from Eq. flMl) that oscillations superimposed upon 
a uniform bouncing universe, as in Eq. fl92|) . lead to additional heating. Similarly, the cases 
where (Af ^) < signal a reduction in quantum uncertainty, or a form of gravitational 
squeezing, analogous to effects near mirrors discussed in Refs. (El, [g], [P]. 

One of the motivations of this study is theoretical, to better understand quantum Brow- 
nian motion in a curved spacetime, as an analog model for the effects of the quantum 
fluctuations of gravity. However, it is also natural to enquire as to whether our results could 
have application to realistic cosmological models. 

One possibility is an additional reheating mechanism after the end of inflation. If inflation 
ends quickly, it is likely that the reheating temperature will exceed the effective temperature 
due to Brownian motion. If reheating is inefficient, however, there is a possibility that 
Brownian motion could play a role. 

Recall that the results in this paper are restricted to the case on non-relativistic motion, 
or when the temperature is small compared to the particle's rest mass energy. This severely 
limit the use of these result for electrons or nucleons. The restriction is less severe for very 
massive particles, such as "wimpzillas" [36|, [13]. These are hypothetical particles with masses 
up to the Planck scale produced at the end of infiation by, for example, gravitational particle 
creation [H, [sH. We plan to extend the study in the present paper to the relativistic motion 
case, which will lift this restriction, and to give a more detailed discussion of applications 
to infiationary cosmology. Another possible extension is to the case of Brownian motion 
produced by fluctuation of non-Abelian gauge flelds. 
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APPENDIX A: MINKOWSKI SPACE CORRELATION FUNCTIONS 

Here we will briefly summarize the calculation of the components of the electromagnetic 
fleld strength tensor correlation function in flat spacetime. Write the Minkowski metric in 
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the form 

ds"^ = -drf + dx^ + dy"^ + dz^ . (Al) 
We are interesting in the components of the field strength tensor correlation function, 

{F^u{x)F^p{x'))m. (A2) 

These are easily computed from the vector potential correlation function using the relation 
Ffj,u = dfj_Ai, — d^Afj,. The vector potential correlation function can, in a suitable gauge, be 
written as 

{A,{x)Mx'))m = ^, (A3) 
where rj^i, is the Minkowski metric tensor, and 

a = i[-(r^ - + (x - x'f + {y - y'f + {z - z'f]. (A4) 
The electric field correlation function, for example, is 

{E^ix) E^{x')) = (F^^x) F-(x')) = - ,.[^l'^tf,2]3 ' (A5) 
where r'^ = — 2Ax'^ and = Aa;^ + Ay^ + Az^. 



APPENDIX B: EVALUATION OF INTEGRALS 



1. Evaluation of Eq. (1321) 

We may evaluate the two integrals in Eq. fl32l) using the residue theorem. First, we 
evaluate the rji integral indicated by: 

-/::-(^) 

In order to do that, let us choose a contour that avoids the third order poles located in 
the real axis at 771 = ±r + 772 (as indicated, in Fig. ([1]), by the letters B and B ). Then, we 
evaluate the integral considering only one of the second order poles, namely, rji = iir^o- In 
Fig. ([T]) we illustrate only the contour in the upper half plane. The contour in the lower half 
plane would give us the same answer, according to the residue theorem. So the first integral 
(/i) has the form 



-{V2 - Vi) - r 
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(Bl) 



J-i — 



Now the second integral is 



CO 



V2 + V0 
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IM A 




t 

FIG. 1: Integration contour for Eq. ()Bip . The lower order poles, which are a pure imaginary 
number, are indicated by the d ot symbol •. The higher order poles has only a real part, as 
indicated by the X symbol and their contribution is a pure imaginary number. 




FIG. 2: This figure illustrated the easiest integration contour we may choose to evaluate Eq. ()B2p . 
It is easiest because we avoid the highest order poles represented by the symbol X in the figure. 
So the contribution of I2 integral is only due to the negative and imaginary lower order pole which 
is indicated by the dot symbol • in the lower half plane. 

The integrand in I2 has a second order pole located at r]2 = ±i?7o and three third and 
one fourth orders poles at 772 = =t^'7o- We evaluate I2 in the lower half plane to avoid the 
highest order poles, as indicated in Fig. ([2]). Even if we had chosen a contour in the upper 
half plane, our answer would be the same, but the way to do that would be harder. So, the 

final result for the mean squared velocity {{Au{ri) Au{r2)) = -^hh) is, 
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{Au{ri)Au{r2)) 



where 6 = — 1. 



r/oe2(12r/2 



r 



+ 



(4r/, 



'2N 



4q^r]oQ 



J2\ 



+ 



(4?7^ + r 



2^4 



2^3 



3r/2e^ 



,/2\ 



,2^2/ 



2^3 



+ 



2(4r/2 + r2)3 
12r/2e 



(4?7^ + r 

{Aril + 
37729(127/2 



2^4 



{Aril + r2)3 



+ 



(477^ + r 
127^20(127/^ 



2\4 



J2\ 



(477, 



„2M 



^27^202 r 72^4Q2 4877^6(127/2 



m2 [ (47/q + r'^Y 



„2^5 



(B3) 



2. Evaluation of Eq. (156]) 

Next we turn to the evaluation of Eq ( !56l) . Define a new Ji integral as 



h 



+00 



drii 



1 



Oq + a\ tanh 



-(^2 - ViY - r"^ 



(B4) 



We see that in the rji integral we have third order poles located in the real axis at: r/i = 
±r + ri2, and an infinite number of poles at tanh(7/i/7/o) = — a^, where = al/af > 1. 
These are the first order poles, and can be written as 



Vik = Vio + 2(2^ + l)^7r?7o, 



(B5) 



with k = 0, ±1, ±2, and 7/10 



^^oln 



Following the contour indicated in Fig. we obtain 
h 



°° 2zr7o cosh 
fc=i 1 



2 / 

Vo 



-{v2 - vikY - 



_7r[-(7/2 - r]ik)^ + '"^l^ 



(B6) 



Note that a contour in the lower half plane would give us the same answer. With this result, 
the 7/2 integral turns into 



k=l 



2ir]o 



cosh^ ( 

^0 



+00 



dr]2 



X 



» Oq + cif tanh 
-(^2 -Vik^ -r"^ 

_7r[-(7/2 -7/ifc)2 + r2]J 



(B7) 
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FIG. 3: One possible contour to integrate Eq. (|B4p . We avoid the higher order poles, indicated 
by X symbol, because their contribution is a pure imaginary number. The contour is, however, 
an infinity contour that encloses the infinity positives lower order poles which are indicated by the 
symbol •. 




FIG. 4: The figure shows the easiest possible contour to integrate Eq. (IB7p . because the infinities 
higher order poles (the X symbol) are excluded. The enclosed poles are the infinities lower orders, 
which are indicated by • symbol. 

It has third order poles, but now located at, 772 = ±r + rjik, and first order poles at 
f]2i = V20 + |(2Z + l)27r?7o, with 7720 = ''7105 and I = 0, ±1, ±2, ...These singularities are located 
in the complex plane as indicated in Fig. (jlj). To avoid the third order poles we choose a 
contour in the lower half plane, however the same answer is required if we had chosen the 
contour in the upper plane. 
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So our result is, 



{Au{ri)Au{r2)) = -4-^ sinh 



2 V - 1 



S, 



where 



k=i i'=i 



s 



^0 



-^^(1,2-^) +^^(2,2 + ^ 



(4-^)^(2,2-;^) + (4 + ^)^(2,2+;^ 



(B8) 



(B9) 



In the second form we transformed the double sum into a single sum by defining j = k + l' — 1 
and / = — /. This sum can be evaluated explicitly to yield 



(BIO) 



where \E'(?t,, x) is the nth Polygamma function (This calculation was done using the algebraic 
program, Maple). 



3. Evaluation of Eq. (ITeH 

Let the Ji integral in the variable rji . We may evaluate it using again the residue theorem. 
The T]i integral has third order poles at rji = ±r + ?72, and single poles at rji = iir^Q. Here, we 
can use again the integration contour in Fig. ([T]). However, t he s ymbol • represents singles 
order poles and not second order poles as in the case of Eq. ( IBll) which corresponds to the 
asymptotically fiat bouncing universe. Thus, the integral Ji is expressed as. 



drji 



-{V2 - Vi f - 



-(^2 - - r'^ 



(Bll) 
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The ri2 integral is given by 



+ 00 



1 



^2^2 + ^2^2 



(B12) 



It has two single poles at r}2 = ±«?7o and two third order poles at 772 = ±r + irjQ. In this 
case we choose a contour in the lower half plane, as in the Fig. ([2]), with the symbol • 
representing now first order poles. This contour avoids the third order poles. Even if we 
choose the conto ur in the upper half plane, our result would be the same. Then, after the 
integrations Eq. ( ]B12|) results in 



{Au{ri)Au{r2)) = ^ 



3?7q — 



(B13) 
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FIG. 5: This figure shows that it is possible to choose one contour where the integral is apparently 
null. 



IM A 




Re 



FIG. 6: This figure shows that Eq. (j93p has at least one imaginary pole located at ryo which can 
be enclosed by the contour of Fig. ([5]) . 

4. Evaluation of Eq. (l93l) 

Next we consider the evaluation of Eq. (1931) . The rji singularities has third order poles 
located at rji = ±r + 772, and first order poles at cos^urji) = —a^rjl — a'^rjQ, where = al/af 

or r/ifc = (a^B{r],k) + ^"'^(^ifc) - l) + iC^k + 1)^, with k = 0,±1,±2, and 

Bijjik = ?7i +?7o)- These singularities are located as indicated in Fig. ([5]), with the X symbol 

now representing the third order poles and the several symbols representing the infinities 

numbers of first order poles. So apparently, this integral is zero because we do not have any 
poles enclosed by one of the possible contours. However, ii rj = ix, we can find two real 
poles in the integrand, because cos^urji) = —cP'J^ — ct'^Vo ^ cosh(u;x) = — a^x^ — a'^rjQ. 

Now we have the picture indicated in Fig. ([0. Then, we have at least one contribution 
due the imaginary first order pole represented there by the symbol •. Thus, the rji integral 
ih) is. 
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FIG. 7: This figure shows the easiest possible contour we can choose to integrate Eq. (jBlSp . Here 
the symbol • represents the first order poles located at ix, the X represents the third order poles, 
and • • •• the infinity first order poles at rjik = In ( oP'B{r]ik) + a'^B{r]ik) — 1) ) + ^{2k + 1)5- 



h = 

and the I2 is defined as: 



2alx — aluj sinh(xc(j) 



IX] 



(B14) 



dri2 



2aQX — a\uj sinh(a;a;) J _^ agr/i + agr/g + o-i cos(ti;?72) 



-3(?72 



IX] 



"1 ^ 



7r[-(r/2 



IX] 
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To treat these poles, we should proceed as in the case of the first integral, as indicated in 
Fig. ([7j). Now the poles are located at 772 = ±r + ix (third order), and at 772 = ±ix (first 
order). Thus, we obtain the result. 



-4 



+ 



2% 



r'^\la^x — alu sinh(Li;x)] 2aQX — afu sinh(u;x) 



-e, 



where G is given by: 
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e 



-{2al{-r' + ix) 



a\uj sinh[ci;(— r' + ix)] }^ 



r{ao(— r + ixY + 00^70 + cos[co'(— r + ix)]}'^ 
2al — afuj"^ cos[co'(— r' + ix)] 
2r{ag(— r + ix)'^ + Ogr/g + af cos[co'(— r + ia;)]}^ 

{2ag(r' + ix) — ajuj sinh[to'(r' + ix)]}^ 
r{al{r + ix)"^ + Ogr/g + a{ cos[co'(r + ia;)]}^ 

— {2ao — a\ijj'^ cos[co'(r' + ix)]} 
2r{ag(r + ix)"^ + agi^g + a\ cos[a;(r + ia;)]}^ 



(B17) 
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In order to check if our answer is correct, let us take Oi ~ and x = rjo. Then we obtain, 

-^2 = 2 4r/l 2 I 213' (B18) 

which is the same answer we found in the bouncing case Eq. (1B13I) . Now, using the smallest 
power terms in Oi and r <^ ?7o, we have the squared coordinate velocity written as, 

and when n/ S> no, (At>(ri) At>(r2)) is given by Eq. (IMI) . 



+ 



3ti;sinh(ti;?7o) + 



cosh(u;r7o 
Vo 



(B19) 
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